Complex solitons with power law behaviour in Bose-Einstein condensates near 

Feshbach resonance 
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Complex, localized stable solitons, characterized by a power law behaviour, are found for a quasi- 
one-dimensional Bose-Einstein condensate near Feshbach resonance. Both dark and bright solitons 
can be excited in the experimentally allowed parameter domain, when two and three-body interac- 
tions are respectively repulsive and attractive. These solutions are obtained for non-zero chemical 
potential, unlike their unstable real counterparts which exist in the limit of vanishing fi. The dark 
solitons travel with constant speed, which is quite different from the Lieb mode, where profiles with 
different speeds, bounded above by sound velocity can exist for specified interaction strengths. 
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The dynamics of non-linear waves in Bose-Einstein 
condensate (BEC) is a subject of immense theoretical and 
experimental interest in current literature. The recent 
observation of dark flj , bright solitons 0, S 0, 01 , soliton 
trains [|| and Faraday waves [7] have given considerable 
impetuous to the investigation of the formation mecha- 
nism and control of various non-linear excitations in the 
quasi- ID scenario Q. The mean field equation, govern- 
ing the dynamics of BEC, is the non-linear Schrodinger 
equation with a harmonic trap. The non-linearity orig- 
inating from the two-body interaction is characterized 
by the s-wave scattering length 'a' which can be con- 
trolled through Feshbach resonance l9| , as also the width 
of the transverse profile For a > 0, elastic interaction 
is repulsive and the BEC is stable. Negative scattering 
length implies attractive interaction, where the conden- 
sate is found to be stable up-to a certain limit of the 
number of atoms l^, U, l2|- 



The three-body interaction can be generally treated as 
a perturbation over the two-body case; it becomes signifi- 
cant for short range and larger scattering length, as is the 
case near Feshbach resonance. For a dense atomic media 
also, three-body interaction plays an important role. It 
is estimated for Rb-BEC that, the real part of the three- 
body interaction term is 10"^ — 10'' times larger than the 
imaginary part 13, 14, 15, T^. Hence, we do not consider 
the three-body recombination here, when the correspond- 
ing coupling constant is imaginary. A number of the- 
oretical studies have been carried out considering three 
body interaction in both three- and quasi-one-dimensions 
[S [S [H, S, EH [12]. Locahzed soliton solutions of 
both elliptic function and power-law type have also been 
investigated 2^ ■ Dark soliton of secant hyper- 
bolic form manifested in purely repulsive three-body in- 
teraction regime relevant for Tonks- Girardeau gas 
[i^ [53, [ill- In this case, the soliton velocity is bounded 
above by sound velocity. Real solitons of both types were 
also analyzed in [i^, where the algebraic one was found 



to exist only in the fi —^ limit and was unstable. 

In this letter, we demonstrate the existence of power- 
law type complex solitons in the presence of repulsive 
two- and attractive three-body interactions. Unlike the 
real case, the obtained dark and bright soliton solutions 
can exist for non-vanishing /x and are stable. The dark 
solitons have a constant velocity determined by the inter- 
action strengths, which is quite different from the Lieb 
mode case j24j. Their profiles can change as a function 
of the parameters of the theory. The corresponding ve- 
locities change from zero to sound velocity. Interestingly, 
in the parameter domain where soliton velocity equals 
sound velocity, it is found that the Bogoliubov dispersion 
is of the quadratic type. For specificity, we consider *^Rb 
with m — 1.44 x 10^^^ Kg and the axial density ctq in the 
range 5.43 x lO^cm^' — 9.67 x lO^cm"-'. The transverse 
trap-frequency is taken as uj± — 2ttx 140 rad/sec and the 
two-body coupling constant 52 = 4.95?i x 10^'^cm'^/sec. 
The three-body interaction coefficient 53 has already 
been estimated 3 IS 31, 13]. The above parame- 
ters allow the present solitons to exist, in the domain of 
53, taking values from — lO^^^cm^/sec to — 10~^^cm^/sec 
(scaled by H). This is in the range of theoretically pre- 
dicted value for ^'^Rb. 

A linear stability analysis using spectral method is car- 
ried out, which shows that the obtained solutions are sta- 
ble against small perturbations in both dark and bright 
soliton regimes. Modulational instability (MI) analysis 
[33, 34, 35J reveals that the parameter regimes relevant 
for the solutions are away from the domain of instability. 

The 3D Gross-Pitaevskii (CP) equation for the wave 
function "^{r^t), with an additional three-body interac- 
tion, is given by 

'^~dt^ -^V^* +{V + .92 1*1' +531*1^ - m) (1) 

where /i is the chemical potential. The cylindrical har- 
monic trap is given hy V ^ mLo\{x'^ + j/^)/2 with a tight 
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transverse confinement. For sufficiently small transverse 
dimension of the cloud, the wave function can be written 



x_+y 
2ai 



) and 



as ■4'ir,t) = f{z,t) (j)o with (/)o = y^^exp(- 

a± = ^/h/(moj±). The longitudinal envelope function 
f{z,t) obeys [lElHS, 



dt 2m dz^ 



+ (ff2|/P+ 531/1' -A^)/, (2) 



where the reduced interaction coefficients are 
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For the space-time independent solution, chemical po- 
tential can be written in terms of the asymptotic density 
ctq; 



M = (52 -f- 53 (To) Co- 



(4) 



real positive quantity. For ^ = 0, one only obtains a real 
soliton [25(1 . Figure [1] shows the density profiles of dark 
and bright solitions for different values of g^. Usually, 
repulsive interaction alone creates dark soliton, whereas 
attractive one results in bright solitons in BEC. As both 
types of forces are present in the present system, one gets 
dark and bright solitons, depending on the values of the 
coupling constants (72 and (73. In Fig. (H (73 is increased 
from dark to bright soliton for a particular value of g2- 
The density profile smoothly transits from dark soliton 
to the bright one. Hence, larger the value of three-body 
interaction, greater is the accumulation of atoms in the 
condensate. Physically it amounts to increasing the local 
density of atoms for going from dark to bright regime. 
This leads to a depletion of atoms in the background. 
The solid line in Fig. [T] corresponds to u = case. Thick 
solid line is the homogeneous background a — ao, where 
u = ±51/520-0/™. 



The superfluid velocity is obtained from the continuity 
equation: 



V = u(l ), 

a 



(5) 



where w = A|| and / = The hydrody- 

namic equation for the density is then, 



2m 



{A^i + 2mu^)a^ +2mMVo.(6) 



A power law ansatz 



aiO = (TO 1 
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1 + De 
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is found to solve Eq. [SI where B and D are given by, 



B = 



3.g2 + 8.g3cro 



2530-0 



D = - 



with u = ±{ 



m (3g2 

5*20-0 



93(Toj 



65"3 

25"3Crg 



(8) 



It is transparent that, non-singular solutions exist only 
when ^3 is negative, i.e., attractive three-body interac- 
tion. The value of g2 should be positive from the reality 
of the soliton velocity, implying repulsive two-body in- 
teraction. As mentioned before, u is a constant for given 
parameter values and density. This situation is quite 
different from the Lieb-mode case, where the soliton ve- 
locity can take different values, bounded above by the 
sound velocity. The obtained solutions can be catego- 
rized into three different classes depending on the values 
of gs for a given (72 and (Tq: (i) A dark solition in the 
range —32/20-0 < 53 < —3g2/8aQ, (ii) a constant back- 
ground for 93 = —3(72/80-0 and (iii) a bright soliton for 
— 3(72/8(To < 53 < — 0.28(72/(To. In these regimes /i is a 
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FIG. 1: The density profiles of soliton solutions for differ- 
ent three-body interactions with 32 = 4.95/i x 10~^^cm'^/sec. 
The obtained dark solitions for gs = —g2/2ao (solid line), 
gs = —0A5g2/cro (dotted line) and bright solitons for g^ = 
—0.32(72 /(TO (small-dashed line), gs — — 0.28(^2 /(to (long- 
dashed line). The thick solid line represents the constant 
background density ao for gs — — 3(/2/8(to. 



The appropriately normalized energy functional 
dfdf* ,92, 
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for the soliton profile. Dark soliton with u = corre- 
sponds to energy E = nhao \J 3(72 (To / (64to) , whereas it 
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goes to zero when the background is uniform. Momen- 
tum of the condensate profile 



i>i 4>2) ■ Here, J is a two-dimensional matrix with Jn 
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' 3(g2 + 2g3cro) 



(10) 



It reaches maximum value (Pmax = Trficro) when = 
—§■2.1 (2(Jo)- Figured] depicts the variation of energy with 
momentum for different three-body interaction strengths. 
Positive momentum is the region of dark soliton, where 
as negative one corresponds to bright soliton. Energy 
and momentum vanish at the transition point a = a^. 
Dispersion graph is stiffer in the bright soliton regimes. 
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FIG. 2: Energy vs momentum for the dark and bright solitons 
for ~g2/2ao < gs < — 0.32(^2 /cro with the same g2 used in 
Fig. [1] Energy and momentum are respectively scaled by 
h^ao^ /m X 10~* and hao- 



The number of atoms in the condensate, normalized to 
vanish at a — ao, 



N = 



dz{aQ — a{z)) 



2fc2„2 \ 1/2 



m\93\ 



93^0 I, 



(11) 



shows that the maximum deficiency of atoms in the dark 
soliton regime is = (67r/i^CTQg2/™)^^^- 

We now analyze the dynamical stability of obtained 
solutions using the spectral method [s^, [s^. A small 
perturbation e'*'*(/)(^) of soliton solution satisfies 



A (y5 = A J 9?, 



(12) 



where is a two-dimensional vector and its components 
are real and imaginary parts of the perturbation: ip — 



J22 = and J12 
operator A are 



- J21 = 1. The elements of the matrix 



An 

A21 
A22 



-52(3 /i + /I) -.93(5 ft + /I + 6 /I) + 

ftu^-2 52/l/2-4<?3/l/2|/P, 

d 

-hu—-2 §2 fif2~4: 53 /i /2 l/P and 



-92UI + 3 /2') -.g3(/i' + 5/2^ + 6 fl fl) + /i, 



where f — {fi + 1/2)- The soliton solution is stable if 
real part of the eigenvalue A is negative. 4>i and 4)2 are 
expanded into a spectral series over 800 modes. This nu- 
merical analysis shows that both bright and dark solitons 
solutions are stable in the entire domain of the solutions. 

It is now worth investigating the issue of modulation 
instability since the three-body interaction is attaractive. 
Phenomenon of modulational instability has been exten- 
sively investigated in literature for BEC 4^, 41, 42|. A 
single component BEC with an attractive atom-atom in- 
teraction, can result in modulational instability, when 
the density of atoms exceeds a certain critical value. We 
assume / = (/o + f)e.xp{i4>), where the infinitesimal fiuc- 
tuation / is given by 

/ = f^cos{Kz - nt) + if2sin{Kz - Qt). (13) 

and K are respectively, the frequency and propa- 
gation constant, of the modulated wave. The above 
transformation produces two sets of equations involv- 
ing fi and /2. Non trivial solutions are obtained only 
if K and satisfy the dispersion relation, 2m^^ = 
K^{h^K^/2m - 4|53|/o^ + 252/0'), where < and 
g~2 > 0. If h^Ky2m < ^3^0^ - 252/0'), it would show 
modulation instability. This condition immediately im- 
plies 52 < 2|53|tTo, which is not in the allowed parameter 
range for the obtained solutions. Thus our solutions are 
modulationally stable. 

In conclusion, complex soliton solutions with power 
law decay have been identified in the quasi-one- 
dimensional GP equation with repulsive two- and at- 
tractive three-body interactions. These solutions, when 
superfluid velocity depends on density, show a slower 
asymptotic decay compared to the elliptic function type 
solutions. We have considered the parameters relevant to 
^^Rb, with the three-body coupling constant in the the- 
oretically predicted range. This opens the possibility of 
observing these complex solitons in realistic BEC. Soli- 
ton velocity is fixed by the strength of the interactions 
and are stable against small perturbations. They are also 
modulationally stable. One would like to study their be- 
haviour in a trap for the purpose of coherent control. 
The analysis of two-soliton sector is also an interesting 
problem, as is the investigation in higher dimensions. 
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